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ABSTRACT 

We present deep Keck spectroscopy for a sample of I%u < 22.5 field early-type galaxies selected 
morphologically in the redshift range 0.56 < z < 1.02 in the Hubble Deep Field North (HDF-N). Using velocity 
dispersions determined from the Keck spectra in conjunction with structural parameters measured from the deep 
WFPC2 images we study the evolution of the M/Lb ratio and the fundamental plane (FP) with redshift. For 
the majority of galaxies the trends observed are very similar to those determined earlier for rich clusters. The 
systematic offset between HDF-N galaxies and cluster galaxies is AlnM/L^ = —0.14 ±0.13, corresponding 
to an age difference of only 16 ± 15% at z = 0.88. However, we find enhanced H5 absorption of equivalent 
width 4.0^0 5 A in the mean spectrum of the ten galaxies, indicating the presence of young stars. We infer that 
the galaxies have composite stellar populations, consisting of a low mass young component in addition to a 
dominating old component. As the bulk of the stellar mass must have formed at z ;> 2 our results argue against 
formation scenarios involving major mergers of gas-rich disk systems at 1 <^z <J 1.5, and we conclude that z — 1 
early-type galaxies were either assembled at higher redshift or in mergers involving little gas. The ubiquitous 
enhanced Balmer lines and the presence of tidal features in two of the galaxies lend some support to the latter 
hypothesis. The main uncertainty in the analysis is the small sample; larger samples are needed to study the 
interplay between the evolution of stellar populations and morphology in detail. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

As continued merging is a central prediction of hierarchical 
assembly models, much effort has been expended in attempt- 
ing to detect a decline with redshift in the abundance of well- 
formed early-type galaxies (e.g., Im et al. 2002, Bell et al. 2003, 
Treu 2003). Present results are somewhat inconclusive, largely 
because of the lack of high resolution imaging data for suffi- 
ciently large samples, the difficulty disentangling luminosity 
evolution and changes in the number density, and the effects 
of large scale structure. The most comprehensive study to date 
is the analysis of a very large color-selected sample by Bell et 
al. (2003), who suggest that the stellar mass in red galaxies may 
have increased by a factor 2-3 between z — 1 and z = 0. 

Two alternative routes to the formation history of field 
spheroidals have emerged in the past few years. Treu et al. 
(1999), van Dokkum et al. (2001) and others have attempted 
to constrain evolution in the fundamental plane (FP) differen- 
tially with respect to that observed in clusters. Differential tests 
are less sensitive to selection effects and provide valuable con- 
straints on environmental trends expected in hierarchical mod- 
els (e.g., Kauffmann 1996). Furthermore, the FP gives informa- 
tion on the mass-to-light (M/L) ratios of galaxies, which can be 
used to interpret the observed evolution of the luminosity func- 
tion in terms of the underlying mass function. Again, current 
results are not yet conclusive. Field and cluster samples with 
z < 0.5 have similar M/L ratios (Kochanek et al. 2000; van 
Dokkum et al. 2001; Treu et al. 2001), but there are indications 
that differences set in at higher redshift (Treu et al. 2002). 

A second development follows the discovery by Menanteau, 
Abraham, & Ellis (2001) that ~25 % of field early-type galax- 
ies in the Hubble Deep Fields (HDFs) show color inhomo- 



geneities such as blue cores, quantified by non-uniformity mea- 
sures across the resolved HST image. Such abnormalities may 
indicate recent (merger-induced) star formation and, if so, their 
study offers a possible route to an in situ mass assembly rate 
(Benson, Ellis, & Menanteau 2002). A key question therefore 
is the extent to which spectroscopic diagnostics of recent ac- 
tivity are seen in distant early-type galaxies (e.g., Barger et al. 
1996). 

In this Letter, we address the spectroscopic properties of a 
representative sample of HST-selected early-type galaxies in 
the HDF-N, with a median redshift z ~0.9. First we investigate 
whether the M/L ratios of field and cluster early-type galax- 
ies diverge beyond z ~ 0.5, as suggested by the work of Treu 
et al. (2002). Secondly, we consider the evidence for recent 
star formation and mergers using spectroscopic diagnostics. As 
the HDF-N represents a small sample, we consider this an ini- 
tial exploration of these issues, demonstrating the value of more 
ambitious surveys now possible with wide field spectrographs 
such as DEIMOS (Faber et al 2002). We assume O m = 0.3 and 
fl\ = 0.7 throughout. 



2. DATA 

2.1. Sample Selection and Spectroscopy 

Early-type galaxies in HDF-N were selected morphologi- 
cally from the catalog of Ellis, Abraham, & Dickinson (2001). 
The reliability of the classifications is discussed in some detail 
by Menanteau et al. (1999). Our present sample comprises a 
subset with 7gi4 < 22.5 and z > 0.5. Spectroscopic redshifts are 
available for almost all of the HDF-N sources to our limit from 
Cohen et al. (2000). 

under 
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Table 1 

Early Type Galaxies in HDF-N Proper 



ID 


HDFN_J 


z 


rtot 
'814 


S/N 


a 


error 




Me 


us 


[O II] 








(Vega) 




(km/s) 




(arcsec) 


(/814) 


(A) 


(A) 


1 


123640.02+621207.4 


1.017 


21.84 


14 


235 


46 


-0.680 


21.79 


5.4 


1.7 


2 


123643.16+621242.2 


0.847 


21.28 


23 


225 


32 


-0.603 


21.44 


1.7 


1.1 


4 


123656.65+621220.2 


0.954 


22.09 


13 


174 


32 


-0.490 


22.91 


2.6 


2.4 


5 


123643.63+621218.3 


0.751 


22.26 


11 










7.0 




7 


123646.14+621246.5 


0.904 


21.89 


12 


244 


48 


-0.635 


22.14 


4.6 


< 1.0 


8 


123650.27+621245.8 


0.680 


20.96 


21 


175 


19 


-0.483 


21.95 


3.3 




13 


123644.38+621133.2 


1.013 


21.06 


20 


286 


29 


-0.220 


23.06 


1.3 


< 1.0 


14 


123700.56+621234.7 


0.562 


20.94 


23 


114 


19 


-0.359 


22.30 


4.7 




16 


123643.81+621142.9 


0.765 


20.43 


18 


161 


27 


-0.423 


21.68 


1.3 


3.0 


17 


123646.35+621404.7 


0.960 


20.84 


26 


182 


29 


-0.897 


20.15 


2.9 


3.9 



Spectroscopic observations were conducted at the Keck I 
telescope using the Low Resolution Imaging Spectrograph 
(Oke et al. 1995) on 2000 March 29-30 and 2001 March 27- 
28. A 7 hour exposure in mediocre seeing was obtained in 
2000 using the 600 lines mm~ 1 grating blazed at 7500 A. Su- 
perior data were obtained in 2001 in good conditions, using 
the 600 lines mm -1 grating blazed at 1 ^m. Between expo- 
sures the galaxies were moved along the slit to facilitate sky 
subtraction. The latter 12 hour exposures provided the bulk 
of the data comprising 10 targets. The instrumental resolution 
Omstr ~ 80km s _1 at 8000 A. The spectral coverage is approx- 
imately 7300-9900 A. The data were reduced following stan- 
dard procedures for dithered multi-slit data; the reduction of 
a very similar dataset is explained in detail in van Dokkum & 
Stanford (2003). 

Internal kinematics were determined by fitting broadened 
template star spectra to the galaxy spectra. Fits were performed 
in real space, since S/N weighting and masking of spectral fea- 
tures is more troublesome in Fourier space (see, e.g., Kelson 
et al. 2000a for details). The template stars cover the wave- 
length range 3200 — 5200 A at high resolution (see van Dokkum 
& Stanford 2003). Velocity dispersions were derived for 9/10 
targets (Table 1). The uncertainties do not include a system- 
atic error of ss 8 % as determined from varying the fitting re- 
gion, template stars, and the continuum filtering. The S/N in 
the spectra ranges from 12 to 26 per A . Studies at low redshift 
have shown that for S/N^ 15 systematic effects begin to domi- 
nate the errors (e.g., J0rgensen, Franx, & Kjaergaard 1995), and 
we find that both the S/N and the formal errors indicate that the 
exposure time was just sufficient for measuring velocity disper- 
sions of these faint galaxies. 

2.2. Structural Parameters 

For the FP analysis, effective radii r e and effective surface 
brightnesses fi e were determined by fitting 2D r 1 / 4 models con- 
volved with the PSF to the 7gi4 images (see, e.g., van Dokkum 
& Franx 1996; Kelson et al. 2000b). The drizzled "Version 
2" data release of the HDF-N (Williams et al. 1996) was used 
for the analysis. We used four well exposed, non-saturated stars 
sampling the full spatial varation of the WFPC2 field to approx- 
imate the true PSFs. The rms range in structural parameters is 
!=s 6 % in r e , and ss 10 % in /j, e . The errors are highly correlated, 
and the rms uncertainty in the product r e /i® 8 which enters the 
FP is « 2 %. Median values of r e and [i e (in units of 7gi4 mag- 
nitudes per square arcsecond) are listed in Table 1 . To compare 



galaxies at different redshifts, surface brightnesses were con- 
verted to a common rest-frame band. For galaxies at z = 0.8 — 1 
the observed 7gi4 band closely corresponds to the rest-frame 
B band; hence we transformed the 7gi4 surface brightnesses to 
rest-frame B using the observed V^m — hu colors, as explained 
in van Dokkum & Franx (1996). 

3. MEAN AGES FROM THE FUNDAMENTAL PLANE 

We determine the mean age of the stellar populations in z — 1 
early-type galaxies from the evolution of the mean M/Lb ratio 
(see Franx 1993). Although derivation of absolute ages from 
morphologically-selected samples is susceptible to progenitor 
bias (van Dokkum & Franx 2001), we can constrain the rela- 
tive ages of field galaxies compared to those of cluster galaxies. 

As demonstrated by, e.g., Treu et al. (2001), the evolving 
M I Lb ratio can be derived from samples spanning a large red- 
shift range by calculating the offsets of individual galaxies 
from the prediction of the locally determined FP. The nearby 
FP of cluster galaxies has the form r e oc a 1 20 I,T 80 in the B 
band (J0rgensen, Franx, & Kjaergaard 1996), which implies 
that M/L oc cr - 50 ^'- 25 oc M 0,25 , assuming that early-type galax- 
ies form a homologous family (Faber et al. 1987). Further as- 
suming that the FP of distant field early-type galaxies has the 
same form as that in nearby clusters, the residual from the FP 
is related to an offset in M/L B ratio through 

AM/L B oca L5( V' 2 V (1) 

(see van Dokkum et al. 2001). 

The AlnM/Lfl ratios of the HDF-N galaxies are shown as a 
function of redshift in Fig. 1(a). Included are cluster galax- 
ies, local field galaxies, and the field sample at 0.15 < z < 
0.55 from van Dokkum et al. (2001). We find AlnM/L B oc 
(— 1.25 ±0.15)z for the full sample of field galaxies, compared 
to AlnM/Ls oc (— 1.12 ± 0.11)z for clusters. Assuming the 
same rate of evolution for both samples and using the biweight 
estimator (Beers, Flynn, & Gebhardt 1990), we find a system- 
atic offset between the HDF-N galaxies and cluster galaxies of 
-0.14 ±0.13 in lnM/L B . Assuming a Salpeter (1955) IMF 
and no systematic difference in metallicity, this offset corre- 
sponds to an age difference of 16 ± 15 % at z = 0.88. The slow 
evolution of field galaxies that we measure strengthens conclu- 
sions from earlier work by Kochanek et al. (2000), Treu et al. 
(2001), and van Dokkum et al. (2001). Treu et al. (2002) find 
stronger evolution of AlnM/L B oc —1. 68^037 for their field 
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Fig. 1. — Panel (a) shows the AlnM/L s ratios of individual early-type galaxies in the HDF-N (solid symbols), relative to the prediction from 
the FP of the nearby Coma cluster. Open circles show field galaxies from Faber et al. (1989) and van Dokkum et al. (2001); dots show cluster 
galaxies from J0rgensen et al. (1996), van Dokkum et al. (1998), and Kelson et al. (2000a). All data except those of Faber et al. were analyzed 
using the same procedures. For reference, the lines indicate stellar formation redshift of 3 (solid) and 1.5 (broken). In panel (b)we show the FP of 
field early-type galaxies edge-on, after correcting the surface brightnesses for the observed luminosity evolution of cluster galaxies. High redshift 
field and cluster galaxies follow very similar trends in both panels. 



sample. This (marginally significant) inconsistency could be 
due to subtle morphological selection effects (see §4), differ- 
ences in the analysis, or other effects. 

In Fig. 1(b) we show the field FP viewed edge-on, with sur- 
face brightnesses corrected according to luminosity evolution 
of cluster galaxies. The FP is remarkably similar to that mea- 
sured in local clusters (solid line). The biweight scatter in the 
HDF-N galaxies is ss 0. 1 in log r e . Larger samples are needed to 
determine the scatter more accurately, and to determine whether 
there is evolution in the tilt. 

The similarity between cluster galaxies and field galaxies is 
remarkable in light of expectations from hierarchical galaxy 
formation models, which predict that field galaxies should be 
substantially younger than cluster galaxies. The Diaferio et al. 
(2001) models predict a difference of AlnM/Ls ~ —0.6 (see 
van Dokkum et al 2001), which is ruled out by the data. 

4. EVIDENCE FOR RECENT MASS GROWTH 

Two galaxies deviate significantly from the general trends in 
Fig. 1, suggesting they have young stellar populations. The 
object with the lowest M /Lb ratio (17) is the only galaxy in 
our sample with a blue core (Menanteau et al. 2001), which is 
likely associated with an active nucleus 5 (Hornschemeier et al. 
2000, Sarajedini et al. 2000). In Fig. 2 we show the ultra-deep 
WFPC2 images of these galaxies, after subtracting the best fit- 
ting r 1 / 4 laws. Remarkably, both deviating objects have mor- 
phological signatures indicative of recent accretion or mergers. 
These features have surface brightness = 24.5 — 25, and 
are detectable only because of the depth of the HDF images. 
None of the other seven galaxies show unambiguous deviations 
from axisymmetry. We conclude that there is tantalizing evi- 
dence for recent mass growth in the population of z — 1 early- 
type galaxies. Note that these anomalous objects were not ex- 
cluded from the field vs. cluster comparison in § 3, although the 

5 Curiously our spectrum shows no evidence for nuclear activity. 



M /Lb ratios of these galaxies are highly uncertain (in particular 
for galaxy 17). 




Fig. 2. — Residual images for two deviant objects in the FP, obtained 
by subtracting the best fitting r 1 ' 4 law profiles (right panels). Both 
show clear non-axisymmetric residual structure, probably of tidal ori- 
gin. The left panels show the pixel-by-pixel color variation. These 
faint structures could only have been reliably detected in the ultra- 
deep HDF data. 



5. EVIDENCE FOR RECENT STAR FORMATION 

Whereas M/Lb ratios are sensitive to the luminosity 
weighted mean age of stellar populations, age-sensitive spec- 
tral diagnostics such as the H8 absorption line provide infor- 
mation on the presence of (small amounts of) young stars (e.g., 
Barger et al. 1996). In Fig. 3 we show the mean spectrum 
of the ten early-type galaxies in our sample, created by nor- 
malizing each spectrum at A lest = 4050 A. The median redshift 
z = 0.90, and the biweight average is zbi = 0.86. The S/N is 
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non-uniform as the full wavelength region is not covered by ev- 
ery galaxy; it is highest (<~ 40) in the region near the Ca II H 
and K break. Shown for comparison is the coadded spectrum 
of 356 red galaxies at 0.30 < z < 0.35 in the Sloan Digital Sky 
Survey (SDSS; Eisenstein et al. 2003). 

The mean spectrum clearly reveals prominent Balmer ab- 
sorption lines in addition to weak [O II] 3727 A emission, in 
striking contrast to the color-selected, lower redshift SDSS 
spectrum. The equivalent widths are Wa(HJ) = 4.0+Q5 an ^ 
Wa([Oii])= — 1.7+i'j, using bandpass definitions of Fisher 
et al. (1998). Uncertainties were derived by bootstrap resam- 
pling; we note that Balmer absorption lines at levels consistent 
with the average value can be detected in almost all individual 
spectra, albeit with lower significance (errors w 1 .0 A; Table 1). 
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Fig. 3. — Coadded rest-frame spectrum for 10 HDF-N early-type 
galaxies at I = 0.90. The Balmer absorption lines are clearly enhanced 
with respect to the (color-selected) z = 0.33 mean SDSS spectrum 
(Eisenstein et al. 2003). The SDSS spectrum is offset by 0.35 in F\. 



6. DISCUSSION 

Our small sample of HDF-N galaxies tells a complex tale: 
the similarity of the M /Lb ratios with those of cluster galaxies 
on one hand contrasts with evidence for recent accretion and 
star formation on the other. The Balmer absorption strength 
of 4.0 A is remarkably high, and exceeds the expectation from 
simple passive evolution of a single age population: taking 
2 < Zform < 3 as the formation epoch of stars in cluster and field 



galaxies (van Dokkum & Franx 2001, and § 3), the Vazdekis 
(1999) models predict Wa(H(5) « 1.0 A for galaxies viewed at 
z ~ 0.9. 

Our results can be explained consistently by postulating that 
a substantial fraction of the HDF-N galaxies experienced some 
recent star formation, but that the bulk of the stellar mass was 
formed at redshifts z J> 2. We used the Vazdekis (1999) mod- 
els to generate composite spectra consisting of a pre-existing, 
2.8 Gyr old population containing the bulk of the stellar mass 
and a low mass young component of varying age. Such models 
naturally produce strong effects in the H8 line and only mod- 
est changes in logM/L B . As an example, a 0.5 Gyr old sec- 
ondary star burst containing 3 % of the total stellar mass gives 
W\(H5) = 3.5 A and AlnM/L B = -0.18, entirely consistent 
with the observations. We note that similar models have been 
proposed for nearby early-type galaxies by Trager et al. (2000). 

It is difficult to unambiguously identify the physical origin 
of recent star formation in the HDF-N galaxies. Tidal features 
detected in two galaxies argue that at least in some cases the 
cause may be a merger or the capture of a satellite. The small 
mass associated with the young component argues against ma- 
jor mergers of gas-rich disk systems in the redshift interval 
1.0 <^ z <^ 1.5. Mergers between bulge-dominated, gas-poor 
systems are much more difficult to constrain, as they would not 
have a large effect on the M/L B ratios and do not develop well 
defined tidal tails (e.g., Barnes 1988). Such mergers have been 
observed in clusters (e.g., van Dokkum et al. 1999), and were 
recently invoked to explain the inferred evolution of the lumi- 
nosity density of red field galaxies (Bell et al. 2003). 

Of course, we cannot yet exclude the possibility that the 
small HDF-N sample is special. The stacked spectrum of 15 
"old" Extremely Red Objects (ERO) presented by Cimatti et al. 
(2002) shows no evidence for enhanced Balmer absorption, 
suggesting significant variation in the properties of z — 1 early- 
type galaxies depending perhaps on how they are selected. Only 
by correlating morphologies, M/L ratios, colors and spectral 
features of large samples is it likely that a quantitative assem- 
bly history for z — 1 early-type galaxies can be determined. 

We thank Felipe Menanteau and Bob Abraham for their sig- 
nificant contributions in the early stages of this project and 
acknowledge discussions with Joshua Barnes, Andrew Ben- 
son, Jarle Brinchmann, Alan Dressier, Jim Gunn and Tommaso 
Treu. We thank the anonymous referee for constructive com- 
ments which improved the manuscript. This work was partly 
supported by STScI grant HST-AR-09541.01-A. The authors 
recognize and acknowledge the cultural role and reverance that 
the summit of Mauna Kea has always had within the indige- 
nous Hawaiian community. We are most fortunate to have the 
opportunity to conduct observations from this mountain. 



REFERENCES 



Barger, A. J., Aragon-Salamanca, A., Ellis, R. S., Couch, W. J., Smail, I., & 

Sharpies, R. M. 1996, MNRAS, 279, 1 
Barnes, J. E. 1988, ApJ, 331, 699 

Beers, T. C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 

Bell, E. F, et al. 2003, ApJ, submitted (astro-ph/0303394) 

Benson, A. J., Ellis, R. S., & Menanteau, F. 2002, MNRAS, 336, 564 

Cimatti, A., et al. 2002, A&A, 381, L68 

Cohen, J. G., et al. 2000, ApJ, 538, 29 

Diaferio, A., Kauffmann, G., Balogh, M. L., White, S. D. M., Schade, D., & 

Ellingson, E. 2001, MNRAS, 323, 999 
Eisenstein, D. J., et al. 2003, ApJ, 585, 694 

Ellis, R. S., Abraham, R. G, & Dickinson, M. 2001, ApJ, 551, 111 



Faber, S. M., Dressier, A., Davies, R. L., Burstein, D., Lynden-Bell, D., 
Terlevich, R., Wegner, G. 1987, Faber, S. M., ed., Nearly Normal Galaxies. 
Springer, New York, p. 175 

Faber, S. M., Wegner, G, Burstein, D., Davies, R. L., Dressier, A., Lynden-Bell, 
D., & Terlevich, R. L. 1989, ApJS, 69, 763 

Fisher, D., Fabricant, D., Franx, M., & van Dokkum, P. 1998, ApJ, 498, 195 

Franx, M. 1993, PASP, 105, 1058 

Hornschemeier, A. E., et al. 2000, ApJ, 541, 49 

Im, M., et al. 2002, ApJ, 571, 136 

J0rgensen, I., Franx, M., & Kjaergaard, P. 1995, MNRAS, 276, 1341 
Kauffmann, G. 1996, MNRAS, 281, 487 



van Dokkum & Ellis 



5 



Kelson, D. D., Illingworth, G. D., van Dokkum, P. G., & Franx, M. 2000a, ApJ, 

531, 159 
— . 2000b, ApJ, 531, 137 
Kochanek, C. S., et al. 2000, ApJ, 543, 131 

Menanteau, F., Abraham, R. G., & Ellis, R. S. 2001, MNRAS, 322, 1 
Menanteau, F, Ellis, R. S., Abraham, R. G, Barger, A. J., & Cowie, L. L. 1999, 

MNRAS, 309, 208 
Oke, J. B., et al. 1995, PASP, 107, 375 
Salpeter, E. E. 1955, ApJ, 121, 161 

Sarajedini, V. L., & Gilliland, R. L., & Phillips, M. M. 2000, AJ, 120, 2825 
Treu, T. 2003, in Clusters of Galaxies: Probes of Cosmological Structure and 

Galaxy Evolution, J. S. Mulchaey, A. Dressier, and A. Oemler, eds., in press 
Treu, T., Stiavelli, M., Berlin, G, Casertano, S., & M0ller, P. 2001, MNRAS, 

326, 237 



Treu, T., Stiavelli, M., Casertano, S., Moller, P., & Bertin, G. 1999, MNRAS, 

308, 1037 
— . 2002, ApJ, 564, L13 

Trager, S. C, Faber, S. M., Worthey, G., & Gonzalez, J. J. 2000, AJ, 120, 165 
van Dokkum, P. G. & Franx, M. 1996, MNRAS, 281, 985 
— . 2001, ApJ, 553, 90 

van Dokkum, P. G, Franx, M., Fabricant, D., Kelson, D. D., & Illingworth, 

G. D. 1999, ApJ, 520, L95 
van Dokkum, P. G, Franx, M., Kelson, D. D., & Illingworth, G. D. 2001, ApJ, 

553, L39 

van Dokkum, P. G. & Stanford, S. A. 2003, ApJ, 585, 78 
Vazdekis, A. 1999, ApJ, 513, 224 
Williams, R. E., et al. 1996, AJ, 112, 1335 



